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ABSTRACT 

The dwarf nova SS Cygni is a close binary star consisting of a K star trans- 
ferring mass to a white dwarf by way of an accretion disk. We have obtained new 
spectroscopic observations of SS Cyg. Fits of synthetic spectra for Roche-lobe- 
filling stars to the absorption-line spectrum of the K star yield the amplitude of 
the K star's radial velocity curve and the mass ratio: Kk = 162.5 ± 1.0 km s -1 
and q = Mx/M w d = 0.685 ± 0.015. The fits also show that the accretion disk 
and white dwarf contribute a fraction / = 0.535 ± 0.075 of the total flux at 
5500 A. Taking the weighted average of our results with previously published 
results obtained using similar techniques, we find (K^-) = 163.7 ± 0.7 km s _1 
and (q) = 0.683 ± 0.012. The orbital light curve of SS Cyg shows an ellipsoidal 
variation diluted by light from the disk and white dwarf. From an analysis of the 
ellipsoidal variations we limit the orbital inclination to the range 45° < i < 56°. 
The derived masses of the K star and white dwarf are Mx = 0.55 ± 0.13 M Q 
and M w d = 0.81 ±0.19 M & , where the uncertainties are dominated by systematic 
errors in the orbital inclination. The K star in SS Cyg is 10% to 50% larger 
than an unevolved star with the same mass and thus does not follow the mass- 
radius relation for Zero-Age Main-Sequence stars; nor does it follow the ZAMS 
mass/ spectral-type relation. Its mass and spectral type are, however, consistent 
with models in which the core hydrogen has been significantly depleted. 

Subject headings: binaries: close, stars: fundamental parameters, individual: SS 
Cygni 
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INTRODUCTION 



As the brightest dwarf nova and one of the brightest cataclysmic variables of any kind, 
SS Cygni has been extensively observed and its properties are thought to be well established. 



Like all cataclysmi c variab 



and a white dwarf (I Warner 



es, it is a binary star consisting of a cool, relatively normal star 
19951) . In deed, it was one of the first cataclysmic variables shown 
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later revised to 0.27513 d by 
of the cool star is K4-5 V f 


Stover et al. 


(1980) and 


Cowlev et a 
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Smith & Dhillon 


1998: 


Beuermann ! 
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The K star fills its 

Roche lobe and is transferring mass through the inner Lagrangian point to an accretion disk 
around the white dwarf. During quiescence the disk has a low viscosity and accumulates 
mass, remaining cool and faint. When the surface density of the disk reaches a critical value, 
the viscosity increases and mass flows rapidly through the disk to the white dwarf, releasing 
gravitational energy. The disk becomes hot and luminous, causing the dwarf nova ou tburst 
( ILasotall200ll ). The outbursts recur roughly every 50 days (jCannizzo &: Matteil Il992h . 



Standard techniques for measuring the masses of binary stars do not work for SS Cyg 
because only the spectrum of the K star is visible and the orbital light curve does not show 
eclipses. Previous measurements of the masses of the stars in SS Cyg have, therefore, been 
forced to rely on dubious assumptions about the properties of the accretion disk or about 
the structure of the K star (see section 5). In this paper we use a technique for measuring 
masses that was first developed for the black hole X-ray binaries and does not require these 
assumptions. The technique has three steps (eg, Casares 2006). First, the amplitude of the 
K star's radial velocity curve Kk and the orbital period P can be combined to give the mass 
function 



M wd sin^ i 



K K P 



(1 + g) 2 2ttG' (1) 

where q = Mk/M W( i is the ratio of the mass of the K star to the mass of the white dwarf and 
% is the orbital inclination. Second, the rotational velocity of the K star can be measured 
from the rotational broadening of its absorption lines. If the rotation period of the K star is 
tidally locked to the orbital period, the ratio of its projected rotational velocity V^sinz to 
Kk is a function of the mass ratio 



V rot sin % 



K 



K 



/(g)- 



(2) 



1 Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the Uni- 
versity of Texas at Austin, the Pennsylvania State University, Stanford University Ludwig-Maximilians- 
Universitat Miinchen, and Georg-August-Universitat Gottingen. 
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Third, the distorted K star causes ellipsoidal variations in the orbital light curve. The 
ellipsoidal variations give another relation between the mass ratio and the orbital inclination. 
Light from the accretion disk and white dwarf can distort the ellipsoidal variations, so it is 
also necessary to measure the amount of contaminating light from these other sources. The 
three relations together yield the dynamical masses of the stars and the dimensions of the 
system. The derived masses depend only weakly on assumed properties of the K star and 
the accretion disk. 

We have obtained new spectroscopic observations of SS Cyg with the Hobby-Eberly 
Telescope (HET). We have measured Kk and q by fitting the spectrograms with model 
spectra calculated using our p rograms for synthesizing t he rotationally broadened spectra of 



cool, Roche-lobe-filling stars (jBitner fc Ro 



bmson 



has been measured by iKjurkchieva et al.l (119991 ) and by IVoloshina fc Khruzinal (120001 ) and 



20061 ). The orbital light curve of SS Cyg 



shows a double-humped ellipsoidal variation diluted by light from the accretion disk. We 
analyze the ellipsoidal variations to place limits on the orbital inclination. Equipped with 
Kk, q, and i, we derive Mk = 0.55 ± 0.13 M & and M w d = 0.81 ± 0.19 Mq . These masses 



are significantly smaller than the masses tabulated in the iRitter &: Kolbl (119981 ) catalog, 
M K = 0.704 M© and M wd = 1.19 M Q . 

Section 2 of this paper describes the spectroscopic observations and section 3 describes 
how we measured Kk and q. In section 4 we describe our analysis of the ellipsoidal variations. 
Section 5 gives a detailed comparison of our results to earlier measurements and discusses 
the reasons for the differences. In the final section we derive the masses of the stars and 
discuss some implications of the lower masses. 



2. OBSERVATIONS 

The spectrograms of SS Cyg were acquired w ith the High Resolution Spectrograph 



(HRS) on the Hobby-Eberly Telescope (jTullI Il998l ). We obtained 23 spectrograms in the 



30 day interval from JD 2452084 to JD 2452114 (June and July 2001), each covering the 
wavelength range from 5300 A to 7000 A at a resolution R = 30,000. The exposure times were 
all 600 seconds, long enough to achieve the necessary signal-to-noise ratio without excessive 
spectral smearing from the changing radial velocity of the secondary star as it progressed 
around its orbit. 

The heliocentric Julian dates at the midpoints of the exposures are given in Table 1. 
Figure [1] shows the eruption light curve of SS Cyg in 2001 from the American Association of 
Variable Star Observers (http:/www. aavso.org/data/lcg). The horizontal bar plotted below 
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the light curve shows when we obtained our s pectrograms. The figur e also shows the eruption 
light curve in 1981 when it was observed by iHessman et al.l (119841 ) . Both sets of data were 
obtained when SS Cyg was in its quiescent state, but in 2001 the quiescent visual magnitude 
was m v = 12.0 — 12.1, about 0.5 magnitudes fainter than in 1981. 

The spectrograms were reduced using standard methods and standard IRAF packages. 
Two instrumental issues are of particular interest for the measurements of SS Cyg. The 
HRS ha s proved to have exce ptionally good radial velocity stability, easily adequate for this 
project (ICochran et al.l 12004] ). On the other hand, the HRS is a fiber-fed echelle spectro- 
graph. Telescope tracking errors, variations in seeing, atmospheric dispersion, and other 
instrumental factors all limit the spectrophotometric accuracy of HRS spectrograms. Nor- 
malizing the continuum required care, especially because the free spectral range of the orders 
is rather small (~70 A). 



As the optica 



1989 



Friend et al. 



near-infrared spectrum of SS Cyg has been published many times (|Echevarria et al 



1990; Martinez-Pais et al. 1994; North et al. 2001, 2002; Webb et al 



2002 



Harrison et al.ll2000l ). we show here only enough of the spectrum to demonstrate the rota- 
tional broadening of its absorption lines. Figure [2] shows the spectrum of SS Cyg near 
6120 A. The spectrum shown is a time-averaged spectrum calculated by shifting the 23 in- 
dividual spectrograms to SS Cyg's systemic velocity and then adding them together. For 
comparison the figure also shows a synthetic spectrum for a single, non-rotating star with 
a K4-5 V spectral type calculated using the ATLAS9 stellar atmosphere program and the 
MOOG spectrum synthesis program. We show in the next section that the extra width of the 
absorption lines in SS Cyg is due to a combination of rotational broadening (~89 km s _1 ), 
orbital smearing during the exposures (< 26 km s _1 ), and instrumental resolution of the 
spectrograph (8.2 km s _1 ). 



3. ANALYSIS OF THE SS CYGNI SPECTRUM 



The Synthetic Spectra: We measured K K , q, and the fraction of the flux contributed by 
the accretion disk and white dwarf by fitting the HET spectroscopy with synthetic spec- 
tra for cool, Roche-lobe-filling stars in close binary systems. We computed the synthetic 
spectra with our LinBrod program, which computes the spectra by summing wavelength- 
depe ndent, velocity-shifted, sp ecific intensities over the distorted surface of the lobe-filling 



star (IBitner &; Robinson! 120061 ). The wavelength-dependent specific intensities are calcu- 



( Sneden 


1973; 


Kurucz 


1993) 
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The synthetic spectra for the K star depend strongly on the orbital period, the orbital 
phase, the mass ratio, and the amplitude of its radial velocity curve; less strongly on the 
star's metallicity and temperature; and only weakly on the orbital inclination, the gravity 
darkening exponent, and other input parameters to ATLAS9 and MOOG. As our primary 
goal was to determine Kk and q, we fixed the orbital period at its known value, we assumed 
the K star has a solar metallicity, and we adopted standard values for the gravity darkening 
exponent and the input parameters to ATLAS9 and MOOG. In preparation for the fits, we 
created large grids of synthetic spectra, two dimensions of the grid covering the possible 
ranges of Kk and q, and a third dimension covering 100 orbital phases. Separate grids 
were calculated for different values of the orbital inclination and effective temperature of the 
lobe-filling star. While we did not fit the temperature and inclination, the separate grids 
were used to determine the sensitivity of Kk and q to those parameters and to measure the 
dilution of the K star spectrum by flux from the accretion disk. To speed the calculation of 
Doppler shifts the synthetic spectra were placed on a logarithmic wavelength scale with a 
spacing of 1 km s -1 . 

The synthetic spectra were convolved with the HRS instrumental profile and with an 
orbital smearing profile. We approximated the instrumental profile with a trapezoid. The 
full width at half maximum (FWHM) of the trapezoid was set equal to 8.2 km s~ x , the 
mean FWHM of a dozen isolated, narrow lines in spectrograms of the Th-Ar wavelength 
calibration lamp. We used a rectangle function to model orbital smearing. The width of the 
rectangle at each orbital phase was determined from AV = K K (2nAt/P) cos27T0, where At 
is the length of the exposure and is the orbital phase. The zero point in phase is the time 
of superior conjunction of the white dwarf. For Kk ~ 162 km s^ 1 and a 10 minute exposure 
time the maximum orbital smearing is ~ 26 km s -1 . 

Determination of K K and q: The synthetic spectra were fitted to the observed spectrograms 
by weighted least squares. As the spectrograms are greatly oversampled in wavelength, we 
were able to assign a weight to each wavelength from the internal scatter of the spectrograms. 
We fitted only those wavelength regions with several strong absorption features. There were 
seven of these regions between 5350 A and 6470 A, each 25 A to 57 A wide and each mostly 
within the central 2/3 of a spectral order. As the regions are narrow, we normalized the 
continua of both the synthetic and the observed spectra to 1.0 by fitting straight lines to the 
continua, a different straight line for each region. 

We fitted all the wavelength regions of all the spectrograms simultaneously, minimizing 
the reduced x 2 °f the fit to find the best values of K K and q, and using the behavior of 
X 2 in the neighborhood of the minimum to estimate the internal errors. There are three 
additional parameters in these fits: the flux contributed by the disk and white dwarf, the 
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syste mic velocity 7, and the orbital phase offset A0, where the A</> is measured with respect 
to the iHessman et al.l Jl984f ] ephemeris: T = HJD 2,444,841.86899 and P = 0.27512973 d. 
The flux from the disk and white dwarf was allowed to run free with a different, independent 
value for each of the seven wavelength regions. As the spectrograms cover the orbit of 
SS Cyg fairly uniformly, we decided to determine 7 and A0 separately from Kk and q. 
Initi al estimates f or 7 a nd A<t> were made f r om fi ts based on initial estimates of Kk and q 
from iBeuermannl (120001 ) and lHessman et al.l (jl984j). Holding these initial estimates for 7 and 
A0 fixed, we then re-fitted the spectra to obtain improved values for K K , q, and the diluting 
flux from the accretion disk. With these improved values for Kk and q, we looped back and 
redetermined 7 and A0. The redetermined values of 7 and A0 did not differ significantly 
from the initial estimates, and so we were done. 

The results are K K = 162.5 ±1.0 km s" 1 , q = 0.685±0.015, A0 = 0.0036 ±0.0010, and 
7 = —16.1 km s . Figure [3] shows the contour plot of \ 2 as a function of Kk and q. The sign 
of A0 is in the s ense that the o b served radial velocities occurred a few minutes later than 
predicted by the IHessman et al.l (11984] ) ephemeris. We do not attach a standard deviation 
to our measurement of 7 because we did not observe radial velocity standards, leaving the 
zero point of our radial velocity syste m uncertain. The 7 velocity is, however, consistent 
with the heliocentric velocity found by iNorth et al.l (120021 ). 7 = —13.1 ± 2.9 km s _1 . 



To investigate the dependence of Kk and q on the effective temperature used to calculate 
the synthetic spectra of the K star, we derived K K and q for two values of T e ff, 4560 K and 
5000 K, the form er corresponding to a K4 V spectral type and the latter to K2 V on the 
Tokunagal (120001 ) temperature scale. We note in passing that the local effective temperature 
varies considerably from point to point over the surface of the star. The temperatures we 
quote are the flux-weighted mean effective temperatures. The dependence on temperature 
is weak. The change of T e ff from 4560 K to 5000 K decreased Kk by only 0.9 km s _1 and 
increased q by 0.005. The values of Kk and q reported throughout this paper correspond to 
4560 K. 

The absorption-line profiles depend on the orbital inclination, so the derived values of 
K K and q also depend on inclination. To investigate the strength of this dependence, we 
re-derived Kk and q for several values of i. The dependence is weak. Changing i from 40° to 
50° changed Kk by less than 0.5 km s _1 and q by about ~ 0.001. Since the absorption-line 
profiles depend on the orbital inc lination, they can be used to determine the inclination, at 
least in principle (IShahbazI Il998l ) . In p ractice the dependence is so weak that fits to the 
spectra do not give reliable inclinations (IBitner fc Robinson! 120061 ). Indeed, for SS Cyg our 
fits would suggest an inclination near or lower than 30°, an impossibly low value since it 
would imply a white dwarf mass much greater than the Chandrasekhar limit. As a result, 
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we do not determine the inclination this way. 



Determination of the Flux from the Accretion Disk and White Dwarf: The effect of extra 
flux from the accretion disk and white dwarf is to dilute the absorption lines in the spectrum 
of the K star. Since the depth of the K star's absorption lines depends on its spectral type, 
the measured amount of flux from the disk and white dwarf depends on the spectral type 
and temperature adopted for the star. The dependence is strong. For T e ff = 4560 K the 
fits yielded / = 0.49, where / is the fraction of the flux at 5500 A coming from the disk 
and white dwarf; and for T e ff = 5000 K the fits yielded / = 0.260. To measure / we must, 
therefore, have a more precise value for the temperature of the K star than we have needed 
up to now. We obtain this temperature in two steps: First we determine the precise spectral 
type of the K star, then we assign a temperature to that spectral type. 



With the exception o f I Joyl (119561 ). who classified the late type star as dG5, and of 
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or K5 V spectral type 


Stover et al. 


1980 




Cowlev et al. 
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Harrison et al. 


2000: 


Webb et al. 



20021 ). The spectral type does not depend on the method nor the wave- 
length region used for the measurement. Neither spectral type is preferred over the ot her. 
Thus, in their reviews of cataclysmic variable secondary stars ISmith &; Dhillonl (119981) as- 
signed a spectral type of K5 V to the secondary of SS Cyg while Beuermannl ( 2000 ) assigned 
a spectral type of K4 V. We adopt K4.5 V with a half subclass uncertainty. T he surface of 



the K star is likely to have regions with much lower temperatures (starspots) (IWebb et al. 
20021 ). Their effect on the spectral type is assumed to be negligible. 



From Table 7.6 in iTokunagal (120001 ) th e temperature corr e spond ing to spectral type 
K4.5 V is 4450 K, while from Table 15.7 in Drilling fc Landoltl (120001 ). the temperature is 
4480 K. Allowing for a half spectral-subclass uncertainty and a possible systematic error of 
100 K due to starspots, heating by irradiation, or an incorrect temperature scale, we take 
T e ff = 4470 ± 140 K as the best estimate of the temperature of the K star. 

As noted above, for T e ff = 4560 K we find / = 0.49, and for T e ff = 5000 K we find 
/ = 0.260. Extrapolating somewhat, we find / = 0.535 ± 0.075 for T eff = 4470 ± 140 K. 
We note that errors in assumed metallicty can affect the line depths in our model LinBrod 
spectra and thus the derived flux contribution from the accretion disk and white dwarf. To 
quantify the magnitude of the effect, we compared results for LinBrod model spectra of 
solar metallicty and [M/H] = -0.5. This is a much larger metallicity error than we expect 
is reasonable. Even so, the derived flux contribution from the white dwarf and accretion 
disk only drops by 20-25% for the lower metallicity spectra. We measured the disk flux in 
each of the seven wavelength regions separately. While the results are rather noisy, the flux 
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contributed by the disk is roughly 20% less at 6400 A than at 5500 A. Other spectroscopic 
measurements of the disk fraction are listed in Table 2 for comparison. With the possible 



exception of the measurement by iNorth et al.l (120021 ) there is surprisingly good agreement 
among them. 

An External Check on the Results: As a final external check on our results one of us (BBB) 
re-reduced the spectrograms of SS Cyg with the Figaro package and re-analyzed them with 
the LINFOR spectral synthesis package (Baschek et al. 1966; Lemke 1997, private com- 
munication), both of which are genetically unrelated to the IRAF and LinBrod packages. 
Using LINFOR we generated a model spectrum for the secondary, including a spherically- 
symmetric rotation profile with v sini = 100 km/s, and cross-correlated this model with each 
observed spectrum to get initial radial velocities. We then adjusted the model rotational ve- 
locity and the abundances of the major line-forming atomic species to minimize the residuals 
between the model and the observations. Finally, we used the model to identify regions of 
continuum and pseudo- continuum, and adjusted the adopted continuum levels for the ob- 
served spectra accordingly. The entire sequence was repeated until the results converged. 
We list the resulting radial velocities in Table 1. The rotational velocity varied from about 
84 to 94 km s _1 , depending on orbital phase, with a mean near V rot smi = 89 km s -1 . This 
method does not yield a reliable mass ratio because of the spherical approximation, and it 
does not yield a measure of the dilution because the dil ution is mimicked by a red uction of 



metallicity; but the radial velocities should be unbiased (IBitner fc Robinson! 120061 ). A fit to 



the velocities yielded = 162.8 km s 1 in excellent agreement with the LinBrod results. 



LIGHT CURVE ANALYSIS AND THE ORBITAL INCLINATION 



The S ource of the Orbital Variations: iKjurkchieva et al.l (119991 ) and IVoloshina &: Khruzina 
( 120001 ) have independently published multicolor mean orbita l light curves for SS C y g and 
the two sets of light curves agree to ~0.01 mag. We use the IVoloshina fc Khruzinal (120001 ) 
data for our analysis because they give the yariances of the individual po ints in the light 
curves. The mean V-band light curve from IVoloshina fc Khruzinal (120001 ). reproduced in 
Figure HI shows a distinctive double-humped modulation with minima near orbital phases 
= 0.0 and 0.5, and unequal maxima near = 0.25 and 0.75, where phase zero is again 
defined to be superior conjunction of the white dwarf. We attribute the double humped 
modulation primarily to ellipsoidal variations of the K star. The best justification for this 
interpretation is the quantitative agreement we find between the observed light curves and 
synthetic light curves based on ellipsoidal variations. The double humps have occasionally 
been attributed to spots on the white dwarf, caused perhaps by accretion onto magnetic 
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poles (IKjurkchieva et al.lll999l ). Since the humps remain at the same orbital phases, this 
interpretation would require the that rotation period of the white dwarf be identical to the 
orbital period. There is no evidence that SS Cyg has a magnetic field strong enough to lock 
on to the K star and enforce synchronism (ie, it is not an AM Her star) and, on the contrary, 
there is evidence that the rotation period of the white dwarf is just a few seconds ( jMauche 
2004 ). Spots on the white dwarf cannot then account for the double humped modulation. 



A pure ellipsoidal variation produces two symmetric humps in a light curve, but in 
SS Cyg the hump at = 0.75 has a larger amplitude th an th e hump at (j) = 0.25. The mul- 



ticolor light curves published IKjurkchieva et al.l (119991 ) and IVoloshina fc Khruzinal (120001 ) 
show that difference between the two humps becomes progressively greater at shorter wave- 
lengths, demonstrating that the extra flux in the larger hump is bluer than the flux in the 
pure ellipsoidal variations. The disks in cataclysmic variables generally have a hot spot on 
their edge where they are hit by a stream of gas transferred from the lobe-filling star, and the 
hot spot oft en produces a single hump in the orbital light curve peaking at phases between 
0.7 and 0.9 (jWarnerlll995l ). Since the extra blue flux in SS Cyg is greatest near phase 0.75, 
we attribute the extra blue flux to a hot spot on the disk. O ur interpretation of the light 
curve is, thus, similar to that of IVoloshina fc Khruzinal (120001 ) . 



A Model for the Light Curve: Our goal is to measure the orbital inclination of SS Cyg 
from the amplitude of the ellipsoidal variations. The ellipsoidal variations are, however, 
diluted by unmodulated light from the accretion disk. This complicates the task because the 
ellipsoidal variations of binaries with low inclinations and low dilution are nearly identical to 
those of binaries with higher inclinations and greater dilution. IVoloshina fc Khruzinal (120001 ) 
attempt to break the degeneracy between inclination and dilution by fitting multicolor light 
curves. For this method to work, the spectral energy distributions of all sources of flux 
must be known. They assumed that all the sources emit black body radiation. This is 
a poor assumption for the optical flux from the qui escent accretion disks of dwarf novae 
(jVrielmann &: Baptistal |2002| ; ISaito fc Baptistal 120061 ) and is not really valid even for the 
flux from the two stars, making the orbital inclination determined from the multicolor fits 
unreliable. 

The degeneracy between the orbital inclination and the disk flux can also be broken if the 
disk flux can be independently measured. The disk flux then becomes an input parameter, 
not a fitted parameter, and we can dispense with theoretical models for the disk flux and 
spectral energy distribution. From the fits to the spectrum described in section 3, we do know 
the V-band flux from the disk, so this is the method we will use to break the degeneracy. 



The computer program we use to calculate light curves for SS Cyg is an updated and 
improved version of the program we have used to calculate the light curves of black hole 
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X-ray binaries (jloannou et al.ll2004j ). Although the program accepts many parameters and 
can synthesize light curves for complicated accretion geometries, just three components are 
needed to fit the light curve of SS Cyg: 



The Ellipsoidal Variations: We use Roche geometry to calculate the shape and local 
effective gravity of the lobe-filling K star. The local effective temperature T e ff (local) 
is calculated from the gravity darkening law T e ff (local) oc T e ff\g\P, where q is the 
local gravity, and (3 is the gravity- darkening exponent. According to I Clare tl (j2000al ) 
(3 ~ 0.10 for main-sequence stars with temperatures near 4500 K. The amplitude 
of the ellipsoidal variations depends only weakly on (3. The angle-averaged fluxes 
in specific filters are calculated b y integrating K urucz model spectra over the filter 
response functions. We adopt the IClaretl (l2000bl ) limb darkening law to convert the 
fluxes to the local specific intensities, which are then integrated over the surface of the 
star to produce the observed brightness. 

Constant Flux from the Disk: The main body of the accretion disk is uneclipsed and 
just adds constant extra flux to the light curve. Any contribution from the white dwarf 
is subsumed in the disk contribution. 

Variable Flux from a Bright Spot on the Disk: The disk has a bright spot painted on 
its outer edge. The spot is visible only when on the side of the accretion disk facing the 
observer and its brightness is modulated by geometric foreshortening. The properties 
of the spot are specified by its position and extent along the disk edge and by its 
brightness. 



The Fits to the Light Curve: Although I Voloshina fc Khruzinal (120001 ) published light curves 
in the U, B, and V bands, we fit only the V-band light curve because the distorting effect 
of the bright spot is least in that filter and because the V band is the only one for which 
our spectroscopy gives a measurement of the diluting flux from the disk. For all fits we use 
Kk = 162 km s _1 and q = 0.69. The deduced orbital inclinations are essentially independent 
of Kk, and changing the mass ratio by Aq = 0.01 changes the deduced orbital inclinations 
by just 0.5°. The errors introduced by uncertainties in Kk and q are, therefore, small 
compared to other sources of error. Since we are fitting just the V-band light curve, the 
deduced inclinations are almost independent of the adopted effective temperature, with just 
a second-order dependence through the limb- darkening and the gravity darkening coefficient. 
For convenience internal to the computer program we used T e ff = 4600 K, not 4470 K. 
Numerical tests showed that the ellipsoidal light curves and the derived inclinations are not 
significantly changed by the difference in temperature. 



- 11 - 



We proceeded by choosing a specific value for the fraction of the V-band flux coming from 
the disk (at <fi = 0.25) and then fitting the synthetic light curve to the data by adjusting 
the remaining parameters to minimize the y 2 of the fit. The four physically-significant 
parameters are the orbital inclination, the location of the spot on the edge of the disk, the 
extent of the spot along the edge, and the total V-band flux from the spot. There are also 
two nuisance parameters: the offset in orbital phase between the synthetic and observed 
light curves, and an overall scale factor. We repeated the fit for many different values of 
the fraction of the flux coming from the disk, producing a table of orbital inclinations and 
standard deviations as a function of the disk flux. 

Figure H] shows a synthetic light curve over-plotted on the observed light curve. The 
synthetic light curve corresponds to an orbital inclination of 52° with 56% of the light coming 
from the accretion disk and a phase offset of A0 = 0.01. The amplitude of the ellipsoidal 
variations - and therefore the deduced orbital inclination - is determined almost entirely by 
the amplitude of the hump at phase 0.25, whereas the properties of the spot are determined 
almost entirely by the hump at phase 0.75. As a result, the orbital inclinations depend only 
weakly on the spot model. Figure [5] shows the deduced orbital inclination as a function of 
the fraction of the V band flux coming from the accretion disk. The best fit inclinations 
are shown by the solid line extending from the lower-left to the upper right in the diagram. 
The dashed lines above and below the solid line show the one standard deviation confidence 
limits for the fits. 

The Orbital Inclination: In the previous section we found that the fraction of the V band flux 
coming from the disk was / = 0.535 ± 0.075. The vertical solid line in Figure [5] corresponds 
to / = 0.535 and the two vertical dashed lines to / = 0.46 and / = 0.61. The allowed values 
of the disk flux and orbital inclination lie in the region enclosed by the four dashed lines 
in the figure. Rounding the artificially sharp corners of the region somewhat, we find that 
the orbital inclination of SS Cyg lies in the interval 45° < % < 56°. The range of possible 
inclinations is rather wide. The width is dominated by the uncertainty in /, which traces 
back to the uncertainty in the temperature of the K star determined from its spectral type; 
other sources of error are minor contributors. 



5. COMPARISON TO PREVIOUS RESULTS 

As the literature on SS Cyg is replete with discrepant measurements of Kk, q, i and the 
masses of its components, and as our measurements lead to yet another set of masses, it is 
appropriate to give a critical review of the literature, comparing our results to the previous 
measurements and analyzing the reasons for the differences. 
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Radial Velocity Curve of the K Star: Table 3 lists all the published measurements 
of Kk- The amplitudes fall into three groups. The first consists of early measurements 
yielding amplitudes near 120 km s _1 . In this group the spectrograms were recorded on 
phot ographic plat e s and velocities were measured using traditional techniques. Beginning 
with IStover et al.l (J1980J) the spectrograms were obtained with digital detectors and the 
velocities were usually measured by cross-correlating or directly fitting a template spectrum 
to the spectrograms. Initially controversial, Stover's measurement was the first of a group 
made between 1980 and 1990 yielding amplitudes near 155 km s _1 . There are two main 
reasons for the difference between these and the earlier results. First, the absorption lines 
from the K star are broadened by rotation and veiled by continuum emission from the 
accretion disk, making line positions difficult to measure on photographic plates. Second, 
heating of the K star by irradiation from the accretion disk during dwarf nova outburst s 
distorts the radial velocity curve, changing the measured amplitude (IRobinson et al.lll986l ). 
While the distortion has not been full y modeled theoretically, the presence of the distortion is 
amply demonstrated by observations (IMartin et al.lll989l ; lEchevarria et al.lll989l ). Observers 
have sinc e become careful to me asure the radial velocity curve of the K star during quiescence, 
although IRobinson et al.l (119861 ) suggested that there might be enough irradiation to bias the 
measured value of Kk significantly even during quiescence. 



The radial velocity amplitudes measured by lMartinez-Pais et al.l (119941 ) and lNorth et al. 



( 120021 ) . 162.5 ± 3 and 165 ± 1 km s _1 respectively, and now 162.5 ±1.0 km s _1 by us, com- 
prise the third groupie While the difference between these amplitudes and the ~ 155 km s _1 
amplitudes of the second group is small, it is statistically significant. We attribute the 



difference to a real change in the SS Cyg system. From ICannizzo fc Matteil (119921 ) the 
mean quiescent visual magnitude of SS Cyg between 1896 and 1992 was m v = 11.8 — 
11.9. The quiescent magnitude fluctuated by a few tenths of a magnitude about the mean 
and appeared t o sho w a ~ 0.1 magnitude secular decrease over the 100 year interval. 
Hessman et al.l (119841 ) observed SS Cyg when its quiescent magnitude was m v = 11.5 — 11.6, 
about 0.3 magnitude brighter than the historical average (see Figured]). AAVSO light curves 
(htt p:/www.aavso.org/data/l cg) show that SS Cyg was at m„ = 11.9 — 12 . whe n observed 
by iMartinez-Pais et al.l (119941 ). at m v w 12.0 when observed by lNorth et al.l (120021 ). and from 
Figured] at m v = 12.0 — 12.1 when observed by us. Succinctly, the recent observations were 
obtained w hen SS Cyg was fainte r than its historical average and muc h fai nter than whe n 
observed by lHessman et al.l (J1984J). Following ICannizzo fc Matteil ( 119921 ) and I Warnerl (119951 ). 
we attribute a lower quiescent brightness to a lower accretion luminosity, which implies lower 



3 km s from Robinson et al 



Martinez-Pais et al.l (Il994l) meas ured 162.5 ± 3 km s 1 from their own data but then adopted 158 ± 

(|l986l ) in the rest of their paper. 
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irradiation of the K star and less distortion of its radial velocity curve. As a result, the values 
of Kr in the third group should be a better representation of the true motion of the K star. 

The Mass Ratio: As the white dwarf in SS Cyg is a minor contributor to the flux 
at optical and infrared wavelengths, its radial velocity curve cannot be directly measured 
and the mass ratio must be deduced by other means. Lacking a better alternative it has 
sometimes been assumed that the wings of the emission lines, which come from high- velocity 
gas near the white dwarf, can serve as a proxy for the white dwarf and track its orbital 
motion. Measurements of the amplitude K em of the emission-line radial velocity curve of 
SS Cyg are shown in Table 3. The weighted average of the amplitudes measured with 
electronic detectors is K em « 92 km s _1 . Taking « K a t, s « 152 — 165 km s" 1 , o n e find s 
q = M K /M wd = 0.56 - 0.60, in accord with the value tabulated bv Imtter fc Kolbl (jl998h . 
Although one of us has previ ously argued th at emission lines can yield good qualitative 
estimates of the mass ratios jRobinsonl fl992} ) . more recent work has demonstrated that 
emission-line velocitie s can yield biased o r simply wrong mass ratios and should not be used 
for accurate work [eg, I Welsh et al.l (119931 )]. 



We have measured the mass ratio from the rotational broadening of the absorption 
lines in the spectrum of the K star. This technique has been widely used to measure the 
masses of b lack holes in b lack hole binary sta rs and, if applied with care, gives reliable 
mass ratios (jCasaresll2006l ). ICowley et al.l (119801 ) used line broadening as a consi s tency check 
on other ways to measure the mass ratio of SS Cyg but iMartinez-Pais et al.l (119941 ) were 
the first to use line broadening as the fundamental way to measure the mass ratio, finding 
q = 0.62 ± 0.05. Their measurement was, however, biased by the adoption of a iGrayl (119921 ) 
rotational broadening profile, which is only appropriate for spherical stars, an d by the use 



of app roximate formulae for converting from rotational velocity to mass ratio. [North et al. 
( 120021 ) also measured the mass ratio from the rotational broadening of the K star's spectrum. 
Their method, outlined in lNorth et al.l (120001 ). is similar to ours and like ours should be free 
of bias. Their result, q = 0.68 ± 0.02, is nearly identical to ours. 

The Orbital Inclination: SS Cyg does not eclipse, so the orbital inclination must be 
measured indirectly. One way has been to assume the K star has some known property 
and then work backwards to the orbital inclination. If one assumes that the K star follows 
the mass/spectral type relation for m ain-sequence stars, its mass is .7 — 0.8 M & and the 
inferred orbital inclination is 38° — 40° (ICowley et al.lll980l ; IWalkerlll98ll ). It is more common 
to assume that the K star follows the main-sequence mass/radius relation. There are many 



1962 



Kiplinger 


1979: 


Stover et al. 


1980: 



45° (IKraft 



19901 ). The properties of lobe-filling stars can, however, depart greatly from main-sequence 
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relations because of nuclear evolution before the onset of mass transfer and the los s of thermal 
equilib rium while mass is being transferred. The theoretical calculations by iKolb et al. 



( 120011 ) show that lobe-filli ng; stars with K4 - K5 spectral types can have masses anywhere in 
the range 0.42 — 0.80 M & . iBaraffe fe Kolbl (120001 ) show that the radii of stars that have used 
up some of the hydrogen in their cores and are transferring mass can be larger than the radii 
of normal m ain-sequence stars with the s ame mass by 50% o r more. The observational data 
collected by lSmith fe Dhillonl (119981 ) and iBeuermannl (120001 ) support these calculations. 



There have been several attempts to estimate the orbital inclination from properties of 
the accretion disk. Co mparing the p r ofile o f the H/3 emission line in SS Cyg to theoretical 
profiles for a thin disk, ICowley et al.l (119801 ) derived an inclination between 30 — 40°. If the 
inner radius of the accretion disk is only slightly greater than the radius of the white dwarf 
and if the white dwarf obeys the usual white dwarf mass/radius relation, one can use the 
inner radius to deduce the mass of the white dwarf and again work backwards to the orbital 
inclination. ICowley et al.l (119801 ) estimated the radius fro m the maximum v e locity extent of 



the wings of the Balmer emission lines, finding i »s 46°. iGiovannelli et al.l (119831 ) assumed 
that the periods of dwarf nova oscillations (DNOs) are equal to the orbital periods of gas at 
the inner edge of the accretion disk. A t that time the DNOs of SS Cyg had been observed to 
range from 7.3 s to 10.9 s, from which IGiovannelli et al.l ( 19831 ) derived i = 40°^2°- Doppler 
tomograms demonstrate, however, that the yelocity fields of accretion disks, includin g the 
accretion disk of SS Cyg, are complicated (IRobinson et al.l Il993l : iNorth et al.l 120011 ). and 
the velocities are likely to be especially complicated near the white dwarf where magnetic 
fields m ay play a role ; and the periods of DNOs are only indirectly related to Keplerian 
periods (jWarnerl 120041 ; iKluzniak et al.l 120051 ) . all of which precludes giving much weight to 
inclinations derived from properties of the accretion disk. 

The analysis of orbital light curves generally yield reliable inclinations. This is clearly 
true for eclipsing binaries but extensive observations of the black hole binaries have shown 
that this is also t rue for orbital inclina tions derived from the ellipsoidal variations of Roche- 
lobe-filling stars (jCharles fe Codl2006l ). Ellipsoidal variations are produced primarily by the 
geometric distortions of the star, which, since the star fills its Roche lobe, are a function only 
of the mass ratio. Large amplitude ellipsoidal variations imply a large inclination, low ampli- 
tude a low inclination. The rapid flickering in the light curve of SS Cyg and its awkwardly- 
long orbital period vitiated earlier attempts to measure its ellipsoidal variations (eg, Honey 
et al. 1989). The variations ar e, however, present and have been observed at both ultra- 



violet and optical wavelengths (|Bartolini et al 



Bruch 1990 



Kjurkchieva et al.lll999l ; IVoloshina fe Khruzina 



1985 



Voloshina 



1986; Lombardi et al. 1987 



20001 ) . The ellipsoidal variations 



of SS Cyg are diluted by flux from the accretion disk and this diluti on must be included in 
the models if fits to the light curves are to give reliable inclinations. IVoloshina fe Khruzina 
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( 120001 ) dealt with this by assuming all components of the SS Cyg system radiate like black 
bodies, finding i = 51° — 56°, but their assumption is questionable, raising concerns about 
the reliability of their result. We have estimated the disk contribution externally, from the 
dilution of the K star's spectrum. Consequently we were able to use a far simpler model 
for the light curve, which leads in turn to a more robust estimate of the orbital inclination. 
Nevertheless, the range of incli nations we find, 45° < % < 56 °, is not incompatible with the 
range of inclinations found by IVoloshina fc Khruzinal (120001 ) . The qualitative reason why 
orbital inclinations derived from the ellipsoidal variations is high is that the observed ellip- 
soidal variations have a large amplitude even though they are heavily diluted by flux from 
the accretion disk. Unless the measured orbital light curves incorrectly represent the true 
light curve, it is difficult to escape the higher inclination. 



DISCUSSION 



The Masses of the Stars: From fits of synthetic spectra to spectrograms of SS Cyg 
obtained with the HET we have found the amplitude of the K star's radial velocity curve 
to be Kk = 162.5 ± 1.0 km s _1 and the ratio of the masses of the two stars to be q = 
Mk/M W( i = 0.685 ± 0.015. From an analysis of the ellipsoidal variations in the orbital light 
curve plus an estimate of the dilution of the ellipsoidal variations by flux from the accretion 
disk, we have limited the orbital inclination to the range 45° < i < 56°. 

Follo wing our discussion in sectio n 5, we see no reason not to accept the measurements 
of Kk by iMartinez-Pais et al.l (119941 ) and iNorth et al.l (120021 ) and the measurement of q 
by iNorth et al.l (120021 ). The best estimates of Kr and q are, then, the weighted average 
of their measurements and ours, (K^) = 163.7 ± 0.7 km s _1 and (q) = 0.683 ± 0.012. 
Adopting these values and taking 45° < i < 56°, we find masses Mk = 0.55 ± 0.13 M and 
M w d = 0.81 ± 0.19 M Q . The uncertainties in the masses are dominated by the uncertainty 
in the orbital inclination, which is itself dominated by systematic errors in the temperature 
of the K star. As a result, it is best to interpret these values as the ranges of possible masses 
and the midpoints of the ranges, not the best masses and their standard deviations. A 
Monte-Carlo error analysis suggests that lower masses within the range are somewhat more 
likely than higher masses. The masses of the K star and the white dwarf cannot be chosen 
independently from within their ranges as they are strongly constrained by the mass ratio 
Mk/M W( i = 0.683 ± 0.012. Kepler's third law now yields the separation of the centers of 
mass of the two stars a = (1.36 ± 0.11) x 10 11 cm, where higher values of a correspond to 
higher masses. 



The Evolutionary State of the K Star: If the mean radius of a Roche lobe is Rl, then 



-16- 



Rl/ci is a function only of the mass ratio, Rl/cl = F(q). Since SS Cyg is a mass transfer 
system, the radius of the K star R k equals the rad ius of its Roche lobe and can be calculated 
from Rx/a = F(q). Adopting the lEggletonl (119831 ) relation for F(q) and using our measure- 
ment of a and M K as a function of i, we find Rk/M k = 1.07 at (i = 45°, M K = 0.68), 
increasing to Rk/Mk = 1-47 at (i = 56°, Mk = 0.42), where the units of R and M are solar 
radius and solar mass. 



From the theoretical models of IChabrier fc Baraffd (119971 ) the mass-radius relation for 
zero age main sequence (ZAMS) stars with masses between 0.4 M and 0.6 M & is R = 
0.95M, or R/M = 0.95, where here also the units are solar radius and solar mass. Thus, 
the K star in SS Cyg is 10% to 50% larger than a single unevolved star with the same 
mass - it does not obey the mass-radius relation for unevolved main-sequence starsjf] Nor 
does the K star follow the ZAMS mass/spectral- type relation for Roche- lobe-filling stars. 
Kolb et al.1 (120011 ) have calculated theoretical evolutionary models for the secondary stars of 
mass-transfer binaries. They find that an unevolved lobe-filling star with spectral type K4-5 
(their "standard sequence" ) should have a mass ~ 0.2 greater than the mass of the K 
star in SS Cyg. Finally, according to lHowell et al.l (120011 ) an unevolved lobe-filling star in a 
binary with SS Cyg's 6.6 hr orbital period should have a mass ~ 0.2 M & greater than that 
of the K star in SS Cyg. 



The Kolb et al.1 (120011 ) models do, however, predict lower masses for K4-5 stars in which 
core hydrogen is depleted but not exhausted. If the mixing length is maintained at a = 1.0, 
the depletion must be extreme to match SS Cyg, X c = 4 x 10~ 4 ; but if the mixing length 
is increased to a = 1.9, the depletion can be more moderate. We interpret our results, 
therefore, as showing that the K star in SS Cyg has significantly depleted its hydrogen. 



Harrison et al.1 (120041 ) showed tha t the CO absorption bands in the infrared spectrum 



of SS Cyg are anomalously weak and lHarrison et al. J2005h found weak CO bands in 9 of 
12 other cataclysmic variables they observed. They attribute the weak CO bands to a low 
C abundance and note that C can be depleted by CNO processing. While the CNO cycle 
can indeed deplete C, all that is really needed is to process C to N; the full CNO cycle is 
not required. In either case C should not be depleted in normal main-sequence stars, so the 
weak CO bands are evidence for nuclear evolution. We take this as supporting evidence that 
the K star in SS Cyg is not an unevolved main-sequence star. 



3 If we were to di scard the orbital inc l inatio n determined from the ellipsoidal variations and insist that 
the K star obeys the Chabrier fc Baraffel ( 1997 ) ZAMS mass-radius relation, we would find Mk = 0.81 M©, 
M w d = 1.18 Mq, and i = 42°. These are stringent upper limits on the masses and lower limit on the orbital 
inclination. 
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The Visual Magnitude of SS Cyg: To calculate the expected V magnitude of SS Cyg, 
we begin with the absolute v isual magnitude of a normal K4.5 V star, which from Table 15.7 
of iDrilling fc Landoltl fepQoh is M v = 7.19. The radius of the K star in SS Cyg is 10% 
to 50% larger than normal. Scaled by its surface area the K star's absolute magnitude 
should be greater than normal by 0.21 to 0.88 magnitude, or My K = 6.64 ± 0.34. The 
trigonometric parallax of SS Cyg is 6.02±0.46 milliarcseconds, leading to a distance modulus 
Am = 6.10 ± 0.16. The apparent visual magnitude of the K star alone should then be 
Vk — 12.74 ± 0.37. This must be augmented by the flux from the accretion disk and white 
dwarf, which we have measured to be a fraction / = 0.535 ± 0.075 of the total V-band flux. 
The extra flux raises the expected V magnitude of the entire system to V — 11.92 ± 0.40, 
which is consistent with the visual magnitude m v = 12.0 — 12.1 measured by the AAVSO at 
the time of our observations. 
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Table 1. Radial Velocities of the Secondary Star in SS Cyg 



Date 


Velocity 


Date 


Velocity 


(HJD - 2452000) 


(km s _1 ) 


(HJD - 2452000) 


(km s _1 ) 


83.81229 


136.2 


104.77608 


104.7 


86.80263 


31.0 


105.76826 


-34.6 


91.80944 


150.7 


106.74746 


-56.2 


94.80816 


108.8 


107.78210 


143.3 


95.79076 


-170.5 


109.76639 


-8.5 


97.76919 


-95.9 


110.77827 


138.2 


98.80362 


-155.5 


111.73070 


-139.7 


99.78387 


144.1 


112.74418 


122.4 


100.75669 


-174.9 


113.72400 


-110.2 


101.78850 


-30.6 


113.76127 


18.6 


102.79008 


121.4 


113.94938 


-174.8 


103.78679 


-168.5 
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Table 2. Measurements of the Contribution of Disk Flux to the 
Total Flux from SS Cygni during Quiescence 



Wavelength 


Fraction of Flux 




Region 


from the Disk 


Reference 


near H/3 


0.50 ±0.03 


Hessman et al. (1984) Quoted 






in Hessman (1988) 


5500 A 


0.56 ±0.03 


Wade (1982) 


5500 A 


0.575 ± 0.075 


This work 


6300 - 6500 A 


0.315 ±0.004 


North et al. (2002) 


6400 - 6700 A 


0.45 


Martinez-Pais et al. (1994) 


7000 - 7600 A 


0.50 


Webb et al. (2002) 
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Table 3. Previous Measurements of the Emission- and 
Absorption-line Radial Velocity Curves of SS Cygni 



K K 


K 




(km s _1 ) 


(km s _1 ) 


Reference 


115 


122 


Jov (1956) 


115 


122 


Walker & Chincarini (1968) 




85 ±4 


KiDlineer (1979) 


120 ±6 


118 ±8 


Cowlev et al. (1980) 


123 ±2 


107 ±2 


Walker (1981) 


153 ±2 


90 ±2 


Stover et al. (1980) 


155 ±2 


96 ±3 


Hessman et al. (1984) 


152 ±2 


96 ±5 


Echevarria et al. (1989) 


155 ±3 


92 ±2 


Friend et al. (1990) 


162 ±3 




Martinez- Pais et al. ( 1994) 


165 ±1 




North et al. (2002) 


162.5 ± 1 




This work 
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Fig. 1. — The eruption light curve of SS Cyg in 1981 and 2001 from the American Association 
of Variable Sta r Observers. Th e hori zontal bar plotted below the light curve in 1981 marks 



the dates when iHessman et al.l (ll984j ) obtained their spectrograms of SS Cyg; the horizontal 
bar in 2001 shows when we obtained our spectrograms. Both sets of data were obtained 
when SS Cyg was in its quiescent state, but SS Cyg was ~0.6 mag fainter during quiescence 
in 2001 than it was in 1984. 

Fig. 2. — The upper panel shows a portion of the average of all our spectrograms of SS Cyg. 
Before the spectrograms were averaged their orbital velocities were removed by shifting the 
spectrograms to SS Cyg's systemic velocity. The lower panel shows a synthetic spectrum for a 
single, non-rotating star with a K4-5 V spectral type for comparison. The extra width of the 
absorption lines in SS Cyg is due to a combination of rotational broadening (~89 km s _1 ), 
orbital smearing during the exposures (< 26 km s _1 ), and instrumental resolution of the 
spectrograph (8.2 km s _1 ). 

Fig. 3. — A contour plot of the y 2 of the fits of the LinBrod synthetic spectra to the 
observed spectrum of SS Cyg as a function of K K and q. The contours of constant x 2 are 
evenly spaced multiples of the minimum x 2 beginning at 1.0001 and separated by 0.0002. 
The best fit values are marked by the cross at q — 0.685 and Kg = 162.5 km s— 1. 

Fig. 4. — T he dots and error bars are m easurements of the V-band orbital light curve of 



SS Cyg from lVoloshina fc Khruzinal (120001 ). Orbital phase = corresponds to spectroscopic 
conjunction with the white dwarf behind the K star. The basic double-humped variation is 
produced by ellipsoidal variations of the K star. The extra amplitude and the asymmetry 
of the hump at = 0.75 is produced by a bright spot on the outer edge of the disk. The 
solid line is a fitted light curve for an orbital inclination i = 52° and 56% of the light at 
= 0.25 coming from the accretion disk. As noted in the text and shown in Figure [5J there 
is a strong positive correlation between the adopted third-light fraction and the best fitting 
orbital inclination. 

Fig. 5. — The correlation between orbital inclination and the fraction / of the V-band flux 
contributed by the accretion disk and white dwarf. The best fitting inclinations for any 
adopted fraction are shown by the solid line extending from the lower-left to the upper right 
in the diagram. The parallel dashed lines show one standard deviation confidence limits 
for the fits. The vertical solid corresponds to / = 0.535, derived from the dilution of the 
absorption lines in the spectrum of the K star; and the two vertical dashed lines correspond 
to / = 0.535 ± 0.075, the confidence limits on f. The derived inclination lies in the interval 
45° < i < 56°. 
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0.66 0.68 0.70 0.72 

Mass Ratio 
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0.3 0.4 0.5 0.6 0.7 

Disk Contribution to V-band Flux 



